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Abstract 


We derive the statistical limit of the spectral autocorrelation function and 
of the survival probability for the indirect photodissociation of molecules in 
the regime of non-overlapping resonances. The results are derived in the 
framework of random matrix theory, and hold more generally for any chaotic 
quantum system that is weakly coupled to the continuum. The “correlation 
hole” that characterizes the spectral autocorrelation in the bound molecule 
diminishes as the typical average total width of a resonance increases. 
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Quantum systems that are classically chaotic are believed to exhibit statistical fluctua¬ 
tions in their spectra and wavefunctions that are universahB These universal properties are 
well-reproduced by the assumption that the Hamiltonians belong to an ensemble of Hamil¬ 
tonians that are consistent with the underlying symmetries, but are otherwise random. Such 
random Hamiltonians are described by random matrix theory (RMT),lii and lead to level 
repulsion, long range correlations in the spectrum, and Porter-Thomas statistics for the 
transition intensities. RMT was initially developed to explain the statistical properties of 
the neutron resonances in the compound nucleus.li However, similar statistical fluctuations 
are observed in a variety of molecular systems. These include fluorescence excitation spectra 
of polyatomic molecules, e.g., the vibrational levels and intensities of nitrogen dioxideSIl and 
pyrazine.i Level repulsion was observed in acetylene,0 and statistical signatures of RMT were 
found in the Stark level-crossing spectra of formaldehyde at the dissociation threshold.! At 
low excitations the molecular Hamiltonian is approximately separable (though anharmonic), 
and the vibrational eigenstates are described by normal modes, allowing the assignment of 
good quantum numbers (“regular” regime). However, at higher excitations the normal 
modes are strongly mixed and the spectrum is better characterized in terms of its statistical 
properties (“chaotic” regime). A transition from regular to chaotic spectrum was observed 
in the vibrational levels of CS 2.0 

Statistical measures of a spectrum include the nearest neighbors level spacing distribution 
and spectral rigidity (e.g. the Mehta-Dyson A 3 statistics). Such measures require accurate 
experimental determination of a complete set of levels.! In practice this is not always possible 
since the finite experimental resolution makes it difficult to determine the statistics of the 
small level spacings. In molecular spectroscopy it is particularly difficult to measure a 
complete set of levels because the density of states is very large at higher excitations. 

A measure of chaos that is less sensitive to experimental resolution is the spectral au¬ 
tocorrelation function and its counterpart in the time domain, the survival probabihty.00 
Since the spectral autocorrelation involves the convolution of the spectrum (i.e. strength 
function) with itself, it is less sensitive to experimental noise. The spectral autocorrelation 
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function is easily constructed from the experimental measurements. For example, using 
lasers one can excite molecules that are vibrationally and rotationally cooled, and measure 
their absorption or fluorescence spectra (i.e., the strength function of the dipole operator). 
The convolution of the measured spectrum with itself is the spectral autocorrelator or the 
power spectrum. The corresponding quantity in the time domain, i.e., the Fourier trans¬ 
form of the power spectrum, is the probability that an initially prepared non-stationary 
state (e.g., the state generated by the operation of the dipole operator on the ground state 
of the molecule) remains in its initial state at a later time. For a chaotic or complex sys¬ 
tem, it was suggested that the power spectrum and the associated survival probability of 
the experimentally prepared state are characterized by the existence of a “correlation hole” 
(i.e. a pronounced local minimum) that originates in the repulsion of energy levels. li A 
correlation hole was observed in the spectrum of highly excited vibrational levels of complex 
molecules like methylglyoxallil and acetylene.lii The correlation hole is a signature of chaos 
and it disappears when the dynamics of the system becomes regular. 

The average survival probability and power spectrum in chaotic systems were evaluated 
in a closed form using RMT for the various classes of Gaussian ensembles,^ and in the 
framework of a scattering modehEi Similar expressions can be obtained by considering dis¬ 
crete time evolution for the circular Dyson ensembles.0 The survival probability (and the 
associated power spectrum) carries information on both spectral and eigenvector statistics.0 
The derivation was restricted to systems that are closed and have discrete spectra. However, 
in many physical situations we are interested in the spectrum at high excitations, where the 
eigenstates are coupled to the continuum and become resonances. The coupling of the quan¬ 
tum system to the continuum is expected to diminish the correlation hole. Experimentally 
one usually observes a quenched correlation hole, an effect attributed to unobserved good 
quantum numbers^l and/or the transition to regular dynamics.00 However, since the cor¬ 
relation hole can also be affected by the coupling of the closed system to its environment, it 
is important to derive the correlation hole quantitatively in this more general situation. 

An important class of processes where coupling to the continuum is important is the 
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photodissociation of molecules. In particular we are interested in indirect photodissociation 
which proceeds through resonances of the excited molecular complex.lil A light pulse excites 
the molecule to a higher electronic state at energies above the dissociation threshold, but 
where a potential barrier hinders the immediate dissociation of the molecule. The barrier is 
usually formed because of the avoided crossing of two diabatic electronic potential surfaces 
(one that is binding and one that is repulsive). Vibrational excitations (below and above 
the barrier) become resonances whose widths depend on the coupling between the binding 
and repulsive manifolds. Such resonances can also be populated in molecular reactions that 
proceed through an excited molecular complex. 

Recently, the universal correlator for the bound-to-continuum strength function was de¬ 
rived for any degree of openness of the quantum system,i§ by using the method of the 
supersymmetric non-linear a modehilii However, the general expression turned out to be 
rather complicated. Here we show that in the regime of isolated resonances (i.e. weak cou¬ 
pling) a much simpler closed form can be derived using random matrix techniques. This 
regime of non-overlapping resonances corresponds to indirect photodissociation at energies 
where the corresponding vibrational resonances he below the barrier in the excited electronic 
potential surface. The width of such resonances is smaller than a typical separation between 
vibrational states by a factor of up to several thousand. The derivation in this paper is 
accomplished in the time domain where the survival probability of the experimentally pre¬ 
pared non-stationary state is averaged over the ensemble. This derivation is analogous to 


that used in Ref. [Tj for closed systems. However, in the case of an open quantum system, 
it is necessary to take into account the decay of the system into the continuum. The prob¬ 
ability for an open chaotic system (that is initially bound) to remain bound at a later time 
was calculated in Ref. The coupling of the system to the continuum is described by the 
average partial widths of the system to decay into each of the open dissociation channels. 
The spectral correlator derived in this paper in the quantum ergodic limit (see Eqs. (pT]) and 
(^) below) depends only on these average partial “dissociation” widths (and the symme¬ 
try class of the Hamiltonian), but is otherwise universal. We hnd that the correlation hole 
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diminishes when the partial widths of the resonances become larger or when the number of 
open channels increases. 

An open system is described by an effective Hamiltonian Tig// = H—ir/2, where H is the 
Hamiltonian of the system when it is uncoupled from the continuum, and F is a matrix de¬ 
scribing the coupling of the system to external open channels.In indirect dissociation, 
H represents that part of the Hamiltonian that includes “binding” potential energy surface 
only, while the open channels describe the possible dissociative states of the molecule. The 
eigenstates | n) of Heff have complex eigenvalues En — iTn/2, describing resonances with 
energies En and widths r„. Bound-to-continuum transitions that proceed through resonance 
states are described by the strength function of the corresponding transition operator 

Tn/2 


a{E) = ^ I (n I T | g) 


( 1 ) 


{E - EnY + Vl/A ■ 

In Eq. d^) each transition intensity (to a state imbedded in the continuum) is weighted by 
a Lorentzian centered at the resonance energy. In indirect photodissociation, the photoab¬ 
sorption spectrum is proportional to Eq. (|^ where T is the dipole operator of the molecule. 
In the following we consider only weakly coupled systems that are characterized by non¬ 
overlapping and narrow resonances (relative to their mean spacing). In this case | n) are 
just the eigenstates of the closed system Hamiltonian H, En are the eigenvalues of H (to 
leading order we ignore shifts in eigenvalues), and the widths are given by 


r = V r 

n / , nc 1 


( 2 ) 


n=l 


where F^g =| y^g p is the partial width of a resonance level n to decay into an open channel 
c (7nc is the partial width amplitude). We assume that there are A open channels. In the 
/^-matrix formulation of reaction theory it is possible to express the partial width amplitude 
as an overlap integral of the eigenstate n and the corresponding channel wavefunction at the 
interface of the internal “interaction” region and the external asymptotic region.@ 

The total strength satisfies the sum rule 


a{E)dE = Y,\{n\T\g)\‘^={g\ T^T \ g) , 


( 3 ) 
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so that the experimentally prepared state T\g) can be normalized by dehning |a) = 
T\g)/^{g \ T^T \ g). The normalized strength function a{E) = cr{E)/ J a{E'dE' can then 
be rewritten as 


^ I (u I a) 

n 


1 rj2 

niE- E^y + Vl/2 • 


(4) 


In the following we shall simply use <j{E) to denote the normalized strength function. 

The spectral autocorrelation function G{uj) is dehned in terms of the strength function 
by0 


/ OO 

a{E)a{E + uj)dE . (5) 

-CXD 

In the time domain we can dehne C{t) as the Fourier transform of (t{E), C{t) = 
jy;^exp{-iEt)a{E)dE. Using (|) we hnd 


C{t) = ^ I (n|a) = {a \ a{t)) , (6) 

n 

where | a(t)) = exp(—iTie//^) | «) is the state at time t that evolves from the experimentally 
prepared non-stationary initial state | a) under the effective Hamiltonian 'Heff- Thus the 
Fourier transform of the spectral autocorrelation G{ijj) is just the survival probability P{t) = 
\G{t)\‘^, i.e. the probability at time t that the system remains at its initially prepared state. 
Notice that for an open system, decay (into the continuum) is included in the time evolution. 
Using (H), the survival probability can be expressed as 

P{t) = ^ I (n I a) 1^ e-rn|t| ^ | (n | a) p| (m | a) p 

n n^m 

We would like to calculate the ensemble average of P{t). In RMT, the eigenvalues and 
eigenvectors are statistically independent. This implies that \ {n \ a) p and F„ (which 
are determined by the eigefunctions) are statistically independent from the spectrum E^- 
Moreover, we assume that the prepared state | a) has no overlap with the channels, so 
that the projection of the eigenstate | n) on | a) and on the channels (i.e. partial width 
amplitudes) are also statistically independent. We conclude that all three quantities |(n|a)p. 
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r„ and En, appearing in Eq. (J^), are statistically independent. Conseqnently, the ensemble 
average on the right hand side of (|^) can be taken over each factor separately. In RMT, the 
eigenvector components are distribnted randomly on the A^-dimensional sphere, and thns 
one hnds for a hxed normalized vector la): 


(n I a) PI {m \ a) P 
I (n I a) P 

A^l (n I a) P 


P 

13 + 2 ’ 
(3 + 2 
N(3 + 2 


( 8 ) 


Here /3 is a constant that depends on the corresponding Ganssian ensemble. For systems that 
conserve time-reversal symmetry the relevant ensemble is the Ganssian orthogonal ensemble 
(GOE) and (3 = 1, while systems that break time-reversal symmetry are characterized by 
the Ganssian nnitary ensemble (GUE) and (3 = 2. Ghaotic nnclear and molecnlar systems 
always have (3=1. The partial width amplitndes can be shown to have a joint Ganssian 
distribntionS 


F(7) oc , 


( 9 ) 


where M^c' = is the channel correlation matrix. Eq. (^) follows from RMT nsing the 
fact that the partial widths can be described as the projection of the internal wavefnnction 
(of the closed system) on the open channels. In the following we assnme, withont loss of 
generality, that all the A channels are nncorrelated. Otherwise, it is always possible to 
transform to the set of eigenchannels. The eigenvalnes of M are jnst the average partial 
widths fc of the eigenchannels. Using T^ = I]c=i I 'Inc P and the known distribntion (H) 
of the partial width amplitndes jnc, it is possible to derive the total width distribntion 
Snch distribntions were observed, for example, in the measured decay rates of 
highly excited vibrational states of formaldehyde in its ground electronic states They 

were also reproduced in numerical simulations of a two-dimensional model chaotic cavity.0 
For the ensemble average of the decay factor in Eq. (|^) we then End 


t 


-fS/2 


3 2 ^ = det fid- 


n 


(i + |n) 


13/2 ■ 


( 10 ) 
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To leading order in 1/iV, r„ and T^ {n 7 ^ m) are statistically independent so that 
e 2 i“i = I e 2 ) . From the known RMT spectral statistics we also have 


{N - l)e^iEn-Em)t = § ^ 


(11) 


JhJ ’ \tH^ 

where tu = 27rp (p is the average level density) is the Heisenberg time. &2,/3('7‘) is the two- 
level form factor, dehned as the Fourier transform of the two-point cluster function T 2 ,/ 3 (s). 
The form factor is known analytically for both the GOE and GUE ensembles! 


52,/3lG = 


[1 - 2|r| + |r| ln(l + 2|r|)] 9{1 - |r|) -F [|r| In (|^) - l] 9{\t\ - 1) {(3 = 1) 


( 12 ) 


(1 - |r|)0(l - |r|) (/3 = 2) 

Gombining (§), (|10|) and ([Tl|) , and measuring time in units of the Heisenberg time (r = 
t/tn), we hnd for the ensemble average of the survival probability 

1 0 + 2^ 1 


iVF(r) = i(T) - M^) n 




(13) 


C (1+ I T I Tjisf /3 y (1 + 2 I r I ' 

where we have used the transmission coefficients Tc. For a weakly coupled system, the 
average partial width Fc is related to the transmission coefficient Tc in channel c by Fc = 
Tc/27rp.ii In the limit of closed system = 0 for all channels and Eq. ([T3|) reduces to Eq. 


(12) of Ref. |4. 

The ensemble average of the spectral autocorrelation G(a;) is found by taking the in¬ 


verse Fourier transform of Eq. 


Dehning the reduced spectral autocorrelation function 


g{u}) = a{E)a{E + (jj) — a{E) /(^{E) = {p/N)G{uj) — a{E) /(t{E) , we hnd the fol¬ 
lowing universal form 


= 2 


0 + 2 


dr cos(27ra;r) J([ ■ 


p/2 


^ ' ' "" (1 + ITct) 

By dehning a function Ep{uj) in the energy domain 

poo 

Ep{u) =2 dr cos{2ttujt) 


dr cos( 27 ra;r) 62/3 ('t) J([ 


(1 + |r,r) 


p/2 


d+RG. 


( 15 ) 


we can rewrite the reduced spectral autocorrelation in the form 































( 16 ) 



There are two interesting limits for Eq. (pTl). In the limit of a closed system all trans¬ 
mission coefficients vanish (Tc = 0), and Eq. (p!4[) reduces to0 



(17) 


The second term in (^) gives rise to the correlation hole discussed in the introduction. 
Another limit for A equivalent open channels (Tc = T for all c) is the limit A —> cx) such 
that AT = K is kept constant. This is equivalent to taking the limit of a large number 
of equivalent channels while keeping the total average width T of the level constant. In 


this limit the integrand in (|T^ becomes e and T^(t) = {k./2) [(vto;)^ -|- {k/2)‘^] ^ is a 
Lorentzian which is independent of /?. Thus we 




( 18 ) 


Eqs. ([I^) and (|T^) provide a closed expression for the universal spectral autocorrelation 
in weakly open chaotic systems, and is the main result of this paper. In Fig. |l|we compare the 


GOE and GUE reduced spectral autocorrelators for an open system coupled to A = 1, 5 and 
10 equivalent channels with a partial transmission coefficient of T = 0.01 in each channel. 
For reference, we show by a dashed line the spectral correlator for a closed system (not 
including the 6 function at ci; = 0). As the number of channels increases, the correlation 
hole shrinks. It is seen that for the same number of open channels, the GOE correlation 
hole is shallower than the GUE hole. Thus the GOE hole disappears faster when opening 
the system to more channels. This could be explained by the stronger level repulsion in the 
GUE case. A similar effect (not shown here) is observed for a fixed number of channels as 
the transmission coefficient increases. Generally, we observe a shallower hole as the total 
mean resonance width increases.^ 

In conclusion, we have derived in closed form the spectral autocorrelation function and 
survival probability for bound-to-continuum transitions that proceed through a regime of 
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isolated resonances. The correlator is universal and depends only on the symmetry class 
(i.e. conserved or broken time-reversal symmetry), and the partial dissociation widths (or 
transmission coefficients) for the open channels. The correlation hole, a signature of chaos 
in both the power spectrum and survival probability, disappears gradually as more channels 
are opened, and it does so faster in systems that conserve time-reversal symmetry. It would 
be interesting to compare our prediction for the correlation hole in open systems (for the case 
of time-reversal symmetry) with actual measurements or numerical simulations of indirect 
photodissociation spectra of complex molecules. 
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FIG. 1. The reduced spectral autocorrelation function g(uj) vs. a; for A = 1, 5 and 10 equivalent 
channels with T = 0.01 (solid lines). The minimum of (“correlation hole”) becomes shallower 
as the number of channels increases. For reference we also show by dashed line g{ijj) for a closed 
system (without the 5 function at w = 0). Top; GOE statistics (/3 = 1); bottom: GUE statistics 
{(3 = 2). 
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